Unified Theoretical Framework for Polycrystalline Pattern Evolution 
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The rate of curvature-driven grain growth in polycrystalline materials is well-known to be limited 
by interface dissipation. We show analytically and by simulations that, for systems forming modu- 
lated phases or non-equilibrium patterns with crystal ordering, growth is limited by bulk dissipation 
associated with lattice translation, which dramatically slows down grain coarsening. We also show 
that bulk dissipation is reduced by thermal noise so that those systems exhibit faster coarsening 
behavior dominated by interface dissipation for high Peierls barrier and high noise. Those results 
provide a unified theoretical framework for understanding and modeling polycrystalline pattern 
evolution in diverse systems over a broad range of length and time scales. 
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Polycrystalline patterns are observed in very diverse 
systems including crystalline solids [1], colloidal systems 
[21 [3] , various spatially modulated phases of macromolec- 
ular systems such as diblock copolymers JUGS, an d non- 
equilibrium (NE) dissipative structures [Bj. When grain 
boundaries (GBs) between domains of different crystal 
orientation are mobile, those patterns generally coarsen 
in time to reduce GB length or area by elimination of 
smaller grains. This coarsening behavior has been exten- 
sively studied because of its practical importance for en- 
gineering polycrystalline materials [7] and its fundamen- 
tal relevance for our general understanding of nonequi- 
librium ordering phenomena. 

The ordering dynamics of modulated phases and NE 
patterns has been extensively studied theoretically [8 - 16 
in the framework of model equations of the form 
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where ip is an order parameter appropriate to each system 
that can be globally conserved (n — 1) or non-conserved 
(n = 0), 77 is a noise uncorrelated in space and time with 
a variance determined by the fluctuation-dissipation re- 
lation (r)(f,t)ri(r',t)) = 2aT(-V 2 ) n 8{r-r , )8{t-t'), and 
J 7 is a Lyapounov functional with a minimum in a lat- 
tice ordered state. Eq. ([!]) has also been proposed as 
a theoretical framework— the phase-field-crystal (PFC) 
model— to study polycrystalline materials on diffusive 
time scales with i/j representing the crystal density field 

PUDS]. While Ec i- has 

been traditionally studied for 
purely relaxational (p = 0) dynamics [5Ml7|. propaga- 
tive (p = 1) wave-like dynamics has also been introduced 
in the PFC framework to mimic phonon-mediated relax- 
ation of the strain field [18 . 

Extensive computational studies of Eq. have shown 
that the characteristic domain or grain size for both 
roll patterns [9HT2] and hexagonal lattices |13fTT6] grows 
~ t q . The exponent q is typically much smaller than 
the q — 1/2 value expected for "normal grain growth" in 



polycrystalline materials [13], and depends on parame- 
ters and noise strength |15| fl6] . While there have been 
theoretical attempts to explain those exponents for roll 
patterns |5HT2"] , the origin of this sluggish (low q) coars- 
ening kinetics is still poorly understood, especially for 
lattices that relate to polycrystalline materials. 

In this letter, we show that the sluggish ordering dy- 
namics of crystal lattices results from the subtle effect 
of "bulk" dissipation linked to lattice translation, absent 
in crystalline solids, but present in Eq. (JT|) that lacks 
Galilean invariance. Internal stresses generated by GB 
motion are known to deform grains and even rotate them 
|21fl2"i] during growth. Our main finding is that the asso- 
ciated dissipation can be, surprisingly, the dominant rate 
limiting mechanism of grain growth for a broad class of 
systems described by Eq. (JT|). This effect is dominant in 
the limit where Peierls-Nabarro (PN) barriers to disloca- 
tion motion are small and curvature-driven growth does 
not require thermal activation. This is the physically rel- 
evant limit for lattice patterns with a wavelength much 
larger than the atomic scale, which coarsen near zero T. 
Interestingly, in the opposite limit of high PN barrier and 
high noise, bulk dissipation plays a subdominant role be- 
cause noise promotes modes of GB motion that reduce 
grain deformation. We summarize here our main results 
and will give a longer exposition elsewhere. 

To highlight the effect of bulk dissipation, we first 
treat analytically the dynamics of a circular grain of ra- 
dius R(t) and misorientation 9(t) with respect to its sur- 
rounding single crystal matrix. A theory for this prob- 
lem has been developed by Cahn and Taylor [32] for 
solids, and basic geometrical aspects of this theory have 
been validated by MD [2J and PFC [23] simulations. 
For small initial misorientation 9(0), grain rotation is 
geometrically necessary under the assumption that the 
number of dislocation along the GB is conserved. Since 
there are rid = 1irR0/b dislocations of Burgers vector 
magnitude 6, this conservation condition implies that 
R(d)6(t) = i?(O)0(O). This rotation can also be inter- 
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Figure 1: Results of isolated grain shrinkage simulations for 
e — 0.12 and T — showing normalized grain area (a) and 
R9/a (b) versus time for different initial grain areas and mis- 
orientations indicated in degrees. Grain rotation is present 
(absent) for misorientation less (larger) than ~ 15°. (c) snap- 
shot of defect field <j> (see text) for the largest area 7° simu- 
lation showing the GB structure, with (11) ((10)) dislocation 
cores appearing as elongated (circular) blue patches, and the 
corresponding ifr field inside the square region (inset). 



preted as a consequence of the geometrical coupling be- 
tween GB motion and a shear stress [22 . 25 , 26] . To treat 
dissipation, we focus for simplicity on non-conversed dy- 
namics (n — 0) but the results hold generally for any n 
or p. For n = and T = 0, Eq. ^implies 



T + E k = - dra{d t xj;) 



(2) 



where the dot denotes total derivative with respect to 
time and = f dr (dtip) 2 /2 is the kinetic energy of 
the ip field. T represents the rate of change of the 



total GB energy of the circular grain, and thus T = 
d[2TrR(t)-f(8(t))]/dt, where 7(0) is the energy per unit 
length of GB. We write the total energy dissipation rate 
as the sum of interface and bulk contributions by split- 
ting the integral on the right-hand-side (r.h.s.) of Eq. 
([2]) as the sum J = fj + f B , where Jj is over a thin an- 
nulus comprising the GB of thickness proportional to the 
dislocation core radius, and f g is over the bulk area of 
the inner grain. Dislocations move radially inward by 
pure glide at a velocity R, and hence \d t ip\ ~ a _1 7? over 
an area a 2 around each dislocation, where a is the lat- 
tice spacing. Therefore, the total interface dissipation 
fjdraidt^) 2 = n d a 2 a(a^ 1 R) 2 /in = a2irRR 2 9/(mb), 
where m is an 0(1) dimensionless pref actor. This yields 
the expression for the mobility M(0) = mb/(a0) for 
-C 1 in agreement with previous studies |2"2l |2"6"] . 
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Figure 2: Comparison of theory (Eq. [4] for the octagonal 
grain with two dislocation types) shown as a black solid (green 
dashed) line for e = 0.12 (e = 0.5), and simulations of Eq. |l]) 
(symbols). Red star symbols are for the modification of Eq. 
|l]) with minimized bulk dissipation (MBD), which causes the 
area decrease rate to become independent of grain size and 
misorientation for small misorientation. 

Next, to compute the bulk dissipation rate 
f B df a(dtip) 2 , we compute the dissipation rate per 
unit area of a crystal field in uniform translation 
at velocity v and then integrate the result over the 
entire grain area. In a region in uniform translation, 
i/j(r,t) fa ipo(f — vt), where ipo(r) is the equilibrium tp- 
field that minimizes and hence dtip = —v ■ Vt/V The 
dissipation rate per area of crystal can therefore be writ- 
ten in the form aKv 2 /A u c , where K = f dr (v ■ Vipo) 2 
is computed over the area A u _ c . of one unit cell (u.c). 
For a square lattice, A u . c , = a 2 and K is independent 
of the direction of v relative to the crystal axes, and 
reduces to K — J® dx'dy' (d x 'ipo) 2 where x' and y' 
are the principal crystal axes. Since v = rO in each 
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small region of a large rotating grain, where r is the 
radial polar coordinates, the total bulk dissipation rate 
is obtained by integrating aKv 2 / 'A u _ c , over the grain 
area: J* 2TrrdraK(r9) 2 / a 2 = airKR^ 2 /(2a 2 ) = 2aE k . 
Combining the above expressions for interface and bulk 
dissipation, Eq. Q becomes 

R<y + R le 9 = -aRR 2 9/{mb) - aKR 4 9 2 /(4a 2 ), (3) 

where we have neglected Ek which can be shown to give 
a negligible higher order contribution (so that the p = 
and p — 1 dynamics in Eq. ([!} yield the same grain ro- 
tation dynamics). Finally, approximating the GB energy 
with a Read-Shockley law 7(6*) = E 9(A C — InO) valid 
for small 9, and using the condition R(t)9(t) = R(Q)9(0) 
that rid is conserved, Eq. ^ yields a single dynamical 
equation for R(t) that can be analytically solved. Its so- 
lution predicts that the grain area (^4 = irR 2 ) decreases 
linearly in time with a rate A = dA/dt given by 

- a 2 /(aA) = Sl a 2 /(mbE ) + s 2 KR(0)9(0)/E , (4) 

where s± — \/(2tt) and S2 — l/(87r), and the first and 
second terms on the r.h.s. correspond to interface and 
bulk dissipation, respectively. Since b ~ a and m and 
K are constants of order unity, Eq. Q predicts that the 
ratio of bulk to interface dissipation is ~ R(0)/a 3> 1, im- 
plying that A is entirely dominated by bulk dissipation, 
which holds for any lattice structure. 




time 

Figure 3: Plot of ordering legnth £s (~ mean grain size) 
vs. time from simulations of Eq. [I] (parameters in legends) 
and minimized bulk dissipation with e — 0.12, a = 8, and 
S = 3.1 10" 3 {S = 6.3 x 10" 3 ) for MBD1 (MBD2); T = 
when not given in legend. Numbers in parentheses are growth 
exponents from straight-line fits at large times (dashed lines). 
Inset: orientation field computed using a wavelet transform 
of ip with a color bar ranging from —45° to 45° [31) . 

To test this prediction, we simulated Eq. (JT|) with p = 
1 and n = and T = j dfuj where |27[ 128] 

u = ip [-e + (V 2 + 1) 2 (V 2 + 2) 2 ] tfj/2 + V//4 - ^, (5) 



has a stable square-lattice ground state that coexists 
with a metastable -0 = liquid-like state for a range of 
H < He- This model was recently shown to produce GBs 
with a similar dislocation content as [001] tilt GBs in 
molecular dynamics (MD) simulations of face-centered- 
cubic bi-crystals [55]. In the analogy with a crystal- 
liquid system, F represents the grand potential that is 
equal in both phases for an equilibrium value jj, = /Ze 
that depends on e [25]. In the simulations reported here 
/i sa 0.9/i£, corresponding to fj, = —0.9 of e = 0.12 and 
/i = —1.67 for e = 0.5. Since the height of the PN bar- 
rier scales ~ exp(— c/e 1 / 2 ) where c is some constant (see, 
e.g. [H]), the height decreases rapidly with decreasing e. 
For e = 0.12, this height is very small so that curvature- 
driven grain growth occurs at T — even for R ~ 100a. 
In contrast, for e = 0.5 the barrier is large and even small 
grains are pinned. For e = 0.5, we carry out simulations 
at finite T and report results for two T values to highlight 
noise effects. Eq. ([lj is simulated using a pseudo-spectral 
method described in [3D] with more details in |31| . 

Fig. [I] shows plots of grain area and R9/a versus time 
together with a snapshot that highlights the structure of 
the GB, consisting of dislocations with Burgers vectors 
described by the Miller indices (11) and (10). Accord- 
ingly, the grain is approximately shaped as an octagon 
with 4 facets made up of [11], [11], [11], and [11] dislo- 
cations, and 4 others with [10], [01], [10], and [01] dislo- 
cations, respectively. Fig. [ip shows that A is constant 
and depends on both initial grain size and misorienta- 
tion, and that R9 is constant. Hence, the number of 
dislocations is conserved until dislocations become too 
closely spaced and annihilate during the final stage of 
grain collapse. In Fig. 2, we compare quantitatively 
A values extracted from linear fits of A vs. time plots 
in the long linear regime before grain collapse to the 
predictions of Eq. Q extended to an octagonal grain 
with two dislocation types: Eq = (Eq 1 + Eq°)/2 and 
(mb)^ 1 = [(toi 6io) _1 + (m 11 6n)~ 1 ]/2 together with the 
shape constants si = l/(4y/2) and S2 = 0.0368 related 
to the perimeter and area of an octagon [31 . This com- 
parison shows an excellent quantitative agreement for 
e = 0.12 and T = for two different a values, confirming 
that inertia (p — 1) is unimportant. The comparison for 
e = 0.5 with large PN barrier shows that bulk dissipation 
is still dominant for intermediate T (T — 0.0325). The 
slope of the curve predicted by Eq. Q fits well the simu- 
lations results, but the curve has a finite intercept at the 
origin corresponding to a finite contribution of interface 
dissipation. This contribution is negligible for e = 0.12 
where the intercept merges with the origin. For e = 0.5 
and larger T (T = 0.065), bulk dissipation is reduced 
and the simulation data is no longer fitted by the theory. 
Analysis of dislocation dynamics in the simulations shows 
that this reduction results from thermally activated dis- 
location reactions. Those reactions reduce the number of 
dislocations , thereby allowing the grain to shrink with 
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less rotation and reducing the contribution of bulk dis- 
sipation relative to interface dissipation. Reduction of 
grain rotation by dislocation reaction is also observed in 
MD simulations of grain rotation |24| . 

To investigate the importance of bulk dissipation in 
grain growth, we simulated Eq. (|T|) with p = 1 and 
n = in large systems of 512 x 512 unit cells with pe- 
riodic boundary conditions. We also simulated a modi- 
fied version of Eq. that minimizes bulk dissipation 
(MBD) by the substitution a — > ah(<fi) where 4>(f) = 
C J df exp(— \f — r*] 2 /2( 2 )\X7ip(r')\ 2 is a smoothly spa- 
tially varying crystalline order parameter (illustrated in 
Fig. JTJd for a rotating grain), with ( = a/2 and the nor- 
malization constant C chosen such that (/> = 1 in ordered 
regions of the lattice and <\> < 1 in disordered regions 
(dislocation cores, GBs, etc). The function h((f) given 
in [31] has limits h(l) = 6 and h(<fi) ~ 0(1) in ordered 
and disordered regions, respectively, thereby allowing us 
to minimize the ratio of bulk to interface dissipation by 
choosing i « 1. Theory for MBD is identical to Eq. 
Q with the second term on the r.h.s. multiplied by S. 
Simulation results in Fig. [2] (red stars) confirm that bulk 
dissipation is negligible for S <C 1. 

We characterized the coarsening behavior by measur- 
ing the ordering scale £5 using the half width 8k at half 
peak of the structure factor S(k,t) = (ip(k,t)ip(—k,t)) 
(fitted to a squared Lorentzian), where the angular brack- 
ets denote averages over all orientations of k in the same 
simulation [101 |3"T]. Plots of = 8k" 1 vs. time in Fig. 
[3] for simulations of Eq. with e = 0.12 and T = 
yield a small coarsening exponent q ss 0.22 as in previous 
studies of hexagonal lattices |13ffT5] . However, MBD sim- 
ulations yield a substantially larger exponent q ss 0.35, 
thereby demonstrating that bulk dissipation associated 
with lattice translation has a significant influence on the 
coarsening kinetics even though grain rotation is more 
constrained in a polycrystal. This conclusion is corrob- 
orated by an analysis of the simulations, which shows 
that grain rotation is more common with MBD and by 
computing the ratio of interface and total dissipation. 
This ratio decreases to a small value with constant a but 
remains approximately unity with MBD [31 . Finally, 
simulations for e = 0.5 show that the growth exponent 
increases markedly with T, as expected from the reduc- 
tion of bulk dissipation by dislocation reactions. More 
important here than the precise values of the exponents, 
which are not exactly constant in time, is the demonstra- 
tion of the key role of bulk dissipation in the ordering 
dynamics. Aspects of the present results, in particular 
the theoretical prediction of the striking dependence of 
grain area shrinkage rate on initial grain size and misori- 
entation (Eq. [4J, should be testable experimentally in a 
wide range of lattice forming systems. 
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